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Scattering  to  extinction  cross-section  ratios,  pse,  were  measured  using  the  NIST  Large  Agglomerate 
Optics  Facility  for  soot  produced  from  ethene  and  acetylene  laminar  diffusion  flames.  Measurements  were 
performed  using  light  sources  at  543.5  nm,  632.8  nm,  and  856  nm.  The  average  scattering  to  extinction 
cross-section  ratios  for  these  wavelengths  are  equal  to  0.245,  0.195,  and  0.195  for  ethene  and  0.311,  0.228, 
and  0.237  for  acetylene.  The  856  nm  measurements  represent  the  longest  wavelength  for  which  accurate 
scattering  measurements  have  been  performed  for  soot.  The  size  distribution  and  fractal  properties  of  the 
two  soots  were  determined  to  assess  the  effects  of  limited  acceptance  angle  range,  finite  size  of  the  sensor, 
and  departure  from  cosine  response  on  the  uncertainty  in  the  measurement  of  pse.  The  expanded  relative 
uncertainty  (95%  confidence  level)  was  found  to  be  ±  6%  at  the  two  visible  wavelengths  and  ±  8%  at  856 
nm.  Both  the  magnitude  and  wavelength  dependence  of  pse  for  the  present  experiments  are  significantly 
different  from  those  reported  by  Krishnan  et  al,  [7]  for  overfire  soot  produced  using  a  turbulent  flame. 
The  results  are  compared  with  the  predictions  of  fractal  optics. 


Introduction 

Two  of  the  defining  characteristics  of  soot  are  its 
absorption  and  scattering  of  electromagnetic  radia¬ 
tion  in  the  visible  wavelength  region.  These  char¬ 
acteristics  affect  the  range  of  atmospheric  visibility 
as  well  as  impair  visibility  in  a  smoke-filled  environ¬ 
ment.  Absorption  and  scattering  properties  of  soot 
are  also  important  in  the  IR  spectrum  for  heat 
transfer  applications,  including  furnaces,  engines, 
and  fires.  There  is  considerable  attention  heeded  in 
the  literature  to  the  measurement  and  prediction  of 
the  extinction  and  absorption  cross-section  of  soot 
[1—4].  There  has  also  been  extensive  study  of  the 
differential  light  scattering  by  soot  in  laboratory, 
scale  fires  [5-7]  as  a  method  for  determining  the 
fractal  dimension  of  the  soot  and  the  size  of  the 
soot  agglomerate  in  terms  of  its  radius  of  gyration, 
Rg.  However,  there  are  few  measurements  of  the 
total  scattering  cross-section  for  flame  generated 
soot. 

Bouguer-Lambert  law  relates  the  light  transmit¬ 
tance,  I/I0,  to  the  mass-specific  extinction  coefficient 
ae  via  the  following  equation: 


—  =  exp(  —  aemsL)  (1) 

lo 

where  ms  is  the  mass  concentration  of  smoke  ag¬ 
glomerates,  and  L  is  the  path  length.  Extinction  is 
the  sum  of  scattering  and  absorption.  Dobbins  et  al. 
[2]  developed  an  approximate  description  of  the  spe¬ 
cific  scattering  and  absorption  of  soot  agglomerates, 
and  the  expressions  [2]  for  these  quantities  are  given 
by 
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where  n1,  n2,  and  Rg  are  the  first  and  second  mo¬ 
ments  of  the  agglomerate  size  probability  distribu¬ 
tion  function,  and  the  mean  of  the  square  of  the 
radius  of  gyration  of  the  soot  agglomerates,  respec¬ 
tively.  The  soot  density  is  ps,  xp  is  nd^/X,  where  1  is 
the  wavelength  of  light,  dp  is  the  soot  primary  par¬ 
ticle  diameter,  D(  is  the  fractal  dimension,  and  the 
functions  of  the  complex  refractive  index  are  defined 
by 
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TURBULENT  COMBUSTION — In-situ  Combustion  Measurements 


Fig.  1.  Schematic  of  the  transmission  cell  reciprocal 
nephelometer,  laminar  burner  system,  and  DOP  aerosol 
generator.  Inset  shows  diagram  of  the  relationship  between 
scattering  pathlength,  dx,  scattering  angle,  6,  and  the  sen¬ 
sor  position  and  orientation. 


trend  suggests  that  scattering  from  soot  maybe  more 
important  than  previously  thought  for  wavelengths 
in  the  near  IR,  at  least  for  large  soot  agglomerates 
typical  of  those  experiencing  long  residence  times 

m. 

Mulholland  and  Choi  [4]  used  a  transmission  cell 
reciprocal  nephelometer  (TCRN),  described  by 
Mulholland  and  Bryner  [3],  to  measure  the  light 
scattering  by  soot  generated  from  laminar  and  tur¬ 
bulent  flames  of  acetylene  and  ethene  at  632.8  nm. 
The  magnitude  of  the  scattering  to  absorption  ratio 
was  approximately  30%  lower  than  that  reported  by 
Krishnan  et  al.  [7]  at  632.8  nm  for  soot  produced 
from  the  same  turbulent  burner  [4], 

The  motivation  for  this  study  is  to  accurately  mea¬ 
sure  pse  and  to  examine  the  influence  of  wavelength 
on  the  light  scattering  by  soot  agglomerates.  In  this 
study,  the  scattering-to-extinction  ratio  at  543.5  nm, 
632.8  nm,  and  856  nm  was  measured  using  the  Na¬ 
tional  Institute  of  Standards  and  Technology  (NIST) 
TCRN  for  soot  produced  from  acetylene  and  ethene 
laminar  diffusion  flames.  Agglomerate  size  distribu¬ 
tions  were  also  obtained  using  a  combination  of 
transmission  electron  microscopy,  optical  micros¬ 
copy  and  image  processing  to  investigate  the  uncer¬ 
tainty  in  the  measurement  of  pse  by  the  TCRN. 


E(m)  =  Im\ 
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Description  of  Experimental  Apparatus 
and  Procedures 


As  pointed  out  by  Dobbins  et  al.  [2],  an  important 
simplification  occurs  for  large  monodisperse  ag¬ 
glomerates  with  4nRg/).  1.  In  this  limit,  the  single 

scattering  albedo  pse ,  defined  as  the  ratio  of  the  scat¬ 
tering  coefficient  to  the  extinction  coefficient,  is 
given  by 

=  cMv/xf~Dl 
^se  6  E(m)  +  Ci(dp/A.)3~Dr 

Cl  =  4F(m)n®-D%(^y  (4) 

For  constant  refractive  index  and  fractal  properties 
with  Df  slightly  less  than  2,  the  value  of  pse  is  pre¬ 
dicted  to  be  roughly  proportional  to  the  primary  par¬ 
ticle  diameter  and  inversely  proportional  to  the 
wavelength. 

Krishnan  et  al.  [7]  measured  pse  for  the  soot  in  the 
overfire  region  of  a  turbulent  burner  for  four  gase¬ 
ous  hydrocarbon  fuels  at  five  different  wavelengths 
in  the  range  between  351  nm  and  632.8  nm.  For  all 
the  fuels,  it  was  found  that  pse  increased  with  wave¬ 
length.  This  result  is  inconsistent  with  the  predic¬ 
tions  of  equation  4  if  the  refractive  index  is  constant. 
In  fact,  Krishnan  et  al.  [7]  attribute  this  increase  to 
a  wavelength  dependence  of  the  complex  refractive 
index,  m  (m  =  n  -I-  ik).  They  point  out  that  this 


Experiments  were  performed  using  the  NIST 
Large  Agglomerate  Optics  Facility  TCRN  to  accu¬ 
rately  measure  optical  properties  of  soot  from  lam¬ 
inar  acetylene  and  ethene  flames.  These  fuels  were 
chosen  because  of  the  large  expected  differences  in 
the  soot  primary  particle  and  agglomerate  sizes.  De¬ 
tailed  descriptions  of  the  apparatus  can  be  found  in 
Mulholland  and  Choi  [4],  Fig.  1  displays  a  schematic 
of  the  experimental  apparatus,  including  the  TCRN 
and  laminar  burner.  The  laminar  burner  fuel  nozzle 
has  an  outer  diameter  (o.d.)  of  12.7  mm,  and  the 
outer  brass  tube  has  an  o.d.  of  10.8  cm.  A  thread  of 
smoke  emitted  by  a  laminar  flame  is  mixed  with  di¬ 
lution  air  as  it  flows  through  a  tripper  plate.  The 
mixture  is  further  diluted  with  air  prior  to  entrance 
into  the  transmission  cell. 

Light  extinction  and  light  scattering  measure¬ 
ments  were  performed  using  a  1  mW  He/Ne  laser 
operating  at  543.5  nm,  a  10  mW  He/Ne  laser  op¬ 
erating  at  632.8  nm,  and  a  30  mW  diode  laser  op¬ 
erating  at  856  nm.  Simultaneous  measurements 
were  performed  using  combinations  of  laser  sources 
(i.e.,  543.5  nm/632.8  nm  and  856  nm/632.8  nm)  in 
order  to  minimize  the  experiment-to-experiment 
variations  in  the  determination  of  scattering  constant 
dependence  on  wavelength.  A  pellicle  beam  splitter 
was  used  (see  Fig.  1)  to  produce  co-finear  beams 


SOOT  SCATTERING  IN  VISIBLE  AND  NEAR-INFRARED 
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Fig.  2.  Laser  transmittance  (solid  line)  and  scattering 
(dashed  line)  measurements  obtained  through  the  DOP 
and  soot  for  laser  wavelengths  of  632.8  nm  and  856  nm. 
Duration  A:  scattered  intensity  and  incident  laser  intensity 
through  clean  air  for  856  nm  laser.  Duration  B:  scattered 
intensity  and  transmitted  laser  intensity  through  DOP  for 
856  nm  laser.  Duration  C:  scattered  intensity  and  trans¬ 
mitted  laser  intensity  through  DOP  for  632.8  nm  laser.  Du¬ 
ration  D:  scattered  intensity  and  incident  laser  intensity 
through  clean  air  for  632.8  nm  laser.  Duration  E:  scattered 
intensity  and  transmitted  laser  intensity  through  soot  for 
856  nm  laser.  Duration  F:  scattered  intensity  and  trans¬ 
mitted  laser  intensity  through  soot  for  632.8  nm  laser. 


through  the  transmission  cell.  The  beams  were  di¬ 
rected  through  the  cell  to  the  silicon  photodiode  de¬ 
tector  using  gold-coated  mirrors.  Rotating  beam 
blocks  were  attached  in  front  of  the  light  sources  to 
selectively  monitor  the  wavelength  of  interest. 

The  scattered  intensity  was  measured  using  a  co¬ 
sine  sensor  and  photo  multiplier  tube  (PMT)  assem¬ 
bly  that  is  positioned  at  the  center  of  the  transmis¬ 
sion  cell  (Fig.  1).  Since  reciprocal  nephelometers 
have  not  been  commonly  used  in  the  combustion 
field,  a  brief  explanation  of  the  principle  of  operation 
is  presented  based  on  the  analysis  of  Gerber  [8,9] 
and  Mulholland  and  Biyuer  [3],  The  inset  in  Fig.  1 
displays  a  schematic  of  the  reciprocal  nephelometer 
with  an  infinitesimal  detector.  The  flux  of  light  scat¬ 
tered  by  the  smoke  particles  from  an  element  of 
length  dx  for  a  laser  beam  with  cross-section  area  as 
in  a  direction  9  toward  the  detector  of  area  is 
given  by  the  following  expression: 

p(xyix  =  (Fpgsad  sm(9)a(9)dx)  ^ 

where  E0  is  the  laser  irradianee,  s  is  the  distance 
between  the  sensor  and  the  element  dx  (where  dx  is 
equal  to  s2d6/D ,  D  is  the  normal  distance  from  de¬ 
tector  surface  to  the  laser  beam),  and  a{0)  is  the 


volumetric  scattering  coefficient  as  a  function  of 
scattering  angle.  The  total  flux  to  the  detector  F,  is 
obtained  by  integrating  over  the  angular  limits  a  and 
/?,  which  approach  0°  and  180°  as  the  length  of  the 
cell  increases. 

F  —  °  s  -  f  a(9)sin(9)d9  (6) 

D  Ja 

The  angular  integral  is  proportional  to  the  orienta¬ 
tion  averaged  scattering  for  non-spherical  particles 
such  as  smoke  agglomerates.  The  above  analysis  is 
for  an  infinitesimal  detector  for  which  eveiy  photon 
incident  on  the  collection  area  is  detected.  For  any 
actual  detector,  there  will  be  specular  reflections  for 
small  angles.  To  have  an  actual  detector  approach 
ideal  behavior,  the  light  must  first  go  through  a  dif¬ 
fuser  with  a  sine  response  relative  to  the  scattering 
angle  or  a  cosine  response  relative  to  the  angle  de¬ 
fined  by  the  normal.  Mulholland  and  Bryner  [3]  find 
from  model  calculations  and  measurements  that 
with  a  small  well-designed  diffuser,  the  TCRN  per¬ 
formance  approaches  that  of  equation  6  given  above. 

In  the  present  experiment  with  D  equal  to  1.3  cm 
and  a  cell  length  of  1  m,  the  collection  angle  of  the 
TCRN  relative  to  the  direction  of  the  laser  beam  was 
from  1.5°  to  178.5°.  The  entire  transmission  cell  was 
covered  with  an  absorbing  black  cloth  to  reduce  am¬ 
bient  and  stray  light  detection.  The  scattering  ex¬ 
periment  was  calibrated  by  simultaneously  measur¬ 
ing  the  light  extinction  and  the  scattering  signal  for 
a  dioctylphthalate  (DOP)  aerosol.  A  constant  flow  of 
aerosol  was  generated  by  nebulizing  a  3%  (by 
weight)  solution  of  DOP  dissolved  in  isopropanol. 
The  aerosol  was  further  diluted  (using  filtered  air) 
prior  to  introduction  into  the  transmission  cell  to 
maintain  a  reasonable  amount  of  attenuation  of  the 
laser  light  within  the  cell.  The  isopropanol  evapo¬ 
rates,  leaving  a  polydisperse  aerosol  with  a  mass 
mean  diameter  of  about  0.5  /nn.  The  key  to  this  cal¬ 
ibration  is  the  fact  that  DOP  does  not  absorb  light 
so  that  the  extinction  constant  measured  for  DOP  is 
equal  to  its  scattering  constant.  The  absolute  scat¬ 
tering  constant  of  soot  was  found  from  the  ratio  of 
the  scattering  signal  measured  for  soot  and  DOP. 
The  incident  light  attenuation  along  the  optical  path- 
length  was  included  in  the  calculation  of  the  scat¬ 
tering  constant. 

The  procedures  for  measuring  the  scattering  con¬ 
stant  for  both  laser  combinations  (543.5  nm/632.8 
nm  and  856  nm/632.8  nm)  are  identical.  Once  a 
steadily  burning  flame  was  established,  clean  air  was 
allowed  to  flow  through  the  transmission  cell  while 
blocking  the  632.8  nm  beam  and  passing  the  856  nm 
beam.  The  scattered  intensity  and  the  incident  laser 
intensity  through  clean  air  were  recorded  during  this 
period  for  the  856  nm  laser  (this  is  denoted  as  A  in 
Fig.  2).  Fig.  2  displays  both  the  transmittance  (solid) 
and  scattered  intensity  (dashed)  measurements. 
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TURBULENT  COMBUSTION — In-situ  Combustion  Measurements 


TABLE  1 

Measurements  of  scattering  to  extinction  ratio  pse  for  soot  produced  from  acetylene  and  ethene  laminar  and 

turbulent  flames 


543.5  nm 

632.8  nm 

856  nm 

ajae  ae  (m2/g) 

^ J  ^ e 

(m2/g) 

ajae  ae  (m2/g) 

Laminar  acetylene,  present  study 

0.311  ±  0.010  8.41 

0.228  ±  0.010 

7.37 

0.237  ±  0.00  5.93 

Laminar  ethene,  present  study 

0.245  ±  0.005  10.18 

0.195  ±  0.003 

8.75 

0.195  ±  0.01  6.28 

Laminar  acetylene11 

0.220  ±  0.007 

7.55 

Laminar  ethenea 

0.186  ±  0.007 

8.50 

Turbulent  acetylene3 

0.250  ±  0.004 

7.80 

Turbulent  ethene8 

0.213  ±  0.010 

8.79 

“Mulholland  and  Choi  [4] . 


Then,  a  constant  flow  of  the  DOP  aerosol  was  intro¬ 
duced  into  the  transmission  cell.  The  resulting  scat¬ 
tered  intensity  and  transmitted  laser  intensity 
through  the  DOP  were  recorded  for  60  s  for  the  856 
nm  laser  (denoted  as  B).  The  856  nm  laser  was  then 
blocked,  and  the  beam  block  in  front  of  the  632.8 
nm  laser  was  removed  to  measure  the  scattered  in¬ 
tensity  and  the  transmitted  laser  intensity  through 
the  DOP  for  the  632.8  nm  laser  (denoted  as  C). 
F resh  air  was  introduced  into  the  cell  to  remove  the 
DOP  aerosol,  and  the  incident  laser  intensity  and 
scattered  intensity  at  632.8  nm  were  measured  (de¬ 
noted  as  D).  Procedure  A  was  repeated  using  the 
856  nm  laser  followed  by  measurement  of  the  scat¬ 
tered  intensity  and  transmitted  laser  intensity 
through  the  soot  (denoted  as  E).  Scattered  intensity 
and  transmitted  laser  intensity  through  the  soot  were 
measured  through  the  soot  using  632.8  nm  laser  (de¬ 
noted  as  F).  Clean  air  was  then  introduced,  and  pro¬ 
cedures  A  through  D  were  repeated  to  confirm  that 
the  magnitude  of  the  drift  in  the  laser  intensities 
during  the  experiment  were  negligible. 

Smoke  samples  were  collected  with  a  thermo- 
phoretic  sampler  for  subsequent  analysis  by  trans¬ 
mission  electron  microscopy  (TEM)  and  optical  mi¬ 
croscopy.  Of  special  interest  was  the  collection  of 
large  agglomerates  to  assess  their  influence  on  the 
total  scattering.  The  large  agglomerates  were  visible 
to  the  unaided  eye  as  distinct  scattering  sources  of 
the  laser  beam.  Smoke  samples  were  obtained  near 
the  center  of  the  transmission  cell.  The  thermo- 
phoretic  sampler  consisted  of  two  brass  plates  sepa¬ 
rated  by  2  mm  which  were  water  cooled  and  could 
be  removed  for  inserting  TEM  grids  or  microscope 
coverslips.  A  nichrome  wire  spaced  0.35  mm  from 
one  of  the  plates  was  maintained  at  100  °C,  and  the 
aerosol  flow  velocity  of  1.2  cm3/s  was  perpendicular 
to  the  direction  of  the  wire.  The  resulting  tempera¬ 
ture  gradient  drives  the  particles  to  the  cooled  sur¬ 
face  where  they  are  collected.  The  band  of  particles 
is  typically  1  mm  wide  and  located  directly  opposite 


of  the  nichrome  wire.  The  sampling  time  of  the  sam¬ 
pler  was  minimized  by  constantly  drawing  a  much 
higher  flow  of  50  cm3/s  from  the  TCRN  with  a  small 
fraction  of  the  flow  going  through  the  precipitator. 


Results 

The  ratio  of  the  scattering  to  extinction  coefficient, 
pse ,  was  calculated  using  both  the  scattering  intensity 
for  DOP  and  soot  and  the  total  extinction  coefficient 
for  DOP  and  soot: 


I^e.dop  — 


HUD 

L 


DOP 


K 


E,soot 


Hip/i) 

L  s 


(7) 


The  PMT  output  voltage  for  scattered  intensity  mea¬ 
surement  for  the  DOP  aerosol  is  Vs  DOP  and  VSsoot 
for  soot.  Since  the  light  extinction  for  DOP  is  com¬ 
prised  entirely  of  the  scattering  component,  the  scat¬ 
tering  coefficient  of  DOP,  Ks  DOP,  is  equal  to  KE  DOP. 
The  attenuation  of  the  laser  beam  through  the  me¬ 
dium  (DOP  or  soot)  will  result  in  a  lower  measured 
light  scattering  coefficient.  Mulholland  and  Bryner 
[3]  obtained  a  correction  factor  for  this  attenuation 
effect  by  comparing  with  a  detailed  numerical  model 
for  the  light  scattering  flux  reaching  the  detector  as 
a  function  of  attenuation  and  particle  size.  The  cor¬ 
rection  factors  for  DOP  and  soot  are  given  as 
exp[KE  DOP  (L/2  +  D)]  and  exp[KE  soot(L/2  +  D)], 
where  D  =  1.3  cm  is  equal  to  the  distance  between 
the  laser  beam  and  the  detector.  The  scattering  co¬ 
efficient  of  soot  was  calculated  from  the  calibration 
with  the  DOP  measurements  as  shown: 


K 


S,soot 


Pse 


VS.soot  exp[KE]SOOt(L/2  +  £>)]  ^ 

Ts.i)0i>('xpA|.;i)0|.  (L/2  +  D)J 

f^S.soot  _  °S 
f^E.soot  °  E 


(8) 


The  average  (based  on  10  measurements  for  each 
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TABLE  2 

Soot  primary  particle  size  and  agglomerate  fractal 
properties 


dp  (nm) 

kf 

Df 

Rg  (nm) 

Ethene 

36.7 

6.85 

1.69 

175 

Acetylene 

50.9 

8.76 

1.61 

241 

10 


0'' 

'ao 
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£0.01 
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0.0001 
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Fig.  3.  Soot  agglomerate  number  distribution  versus 
geometric  mean  radius  of  gyration  for  soot  generated  by 
the  burning  of  acetylene  and  ethene  in  a  laminar  diffusion 
flame. 


fuel  and  wavelength  combination)  scattering  to  ex¬ 
tinction  ratio,  psa,  and  the  specific  extinction  coeffi¬ 
cients  are  summarized  in  Table  1. 

Soot  agglomerate  structure,  including  the  size  dis¬ 
tribution  of  the  primary  particles,  the  fractal  dimen¬ 
sion,  the  number  of  primary  spheres,  and  the  radius 
of  gyration,  was  analyzed  using  TEM.  These  quan¬ 
tities  are  related  by  the  fractal  equation  for  soot  ag¬ 
glomerates  given  by 


An  image  processing  algorithm  was  used  to  measure 
n,  d„,  Rg,  D[,  and  kf  from  digitized  images  of  the  soot 
agglomerate  micrographs  [10],  Measured  values  for 
ethene  and  acetylene  soot  are  shown  in  Table  2. 

A  major  challenge  in  any  total  scattering  instru¬ 
ment  is  the  amount  of  scattered  light  outside  of  the 
acceptance  angle  range  of  the  sensor.  We  have  fo¬ 
cused  on  measuring  the  large  agglomerate  portion 
of  the  size  distribution  since  a  greater  fraction  of  the 
light  scattered  by  larger  agglomerates  is  in  the  for¬ 
ward  direction  and  wall  affect  the  detector  perfor¬ 
mance.  The  overall  distribution  analysis  required  a 


combination  of  TEM  and  optical  microscopy.  The 
lower  magnification  and  resulting  larger  field  of  view 
of  optical  microscopy  are  needed  for  obtaining  a  sta¬ 
tistically  meaningful  number  of  the  infrequent  large 
agglomerates  with  lengths  greater  than  about  3  pm. 
A  limitation  of  TEM  analysis  for  larger  agglomerates 
is  the  enhanced  deposition  of  agglomerates  larger 
than  about  5  /mi  on  the  grid  bars  of  the  TEM  sam¬ 
ples.  In  the  combined  TEM  and  optical  microscopy 
method,  size  classes  of  agglomerates  are  divided 
with  equal  spacing  in  terms  of  the  logarithm  of  ag¬ 
glomerate  size  as  is  widely  done  for  broad  aerosol 
size  distributions  [11].  The  range  in  size  class  is 
maintained  at  a  factor  of  2;  for  example,  all  the  ag¬ 
glomerates  in  the  size  class  between  1  pm  and  2  pm 
are  in  one  class,  those  2  pm  to  4  /an  are  in  another 
class,  et  cetera.  Data  were  obtained  at  a  magnifica¬ 
tion  of  X  2000  with  a  TEM  and  at  X  1000  and  X  500 
with  an  optical  microscope.  The  agglomerates  were 
sized  in  terms  of  (LW)1/2  where  L  and  W  are  the 
length  and  width  of  the  agglomerate.  At  least  20  ag¬ 
glomerates  were  sized  in  each  size  class  to  assure 
reliable  statistics.  Approximately  1500  agglomerates 
were  sized  by  the  TEM  and  another  1500  by  optical 
microscopy  for  both  the  acetylene  soot  and  the  eth¬ 
ene  soot.  The  optical  microscope  data  are  scaled  so 
that  the  TEM  and  optical  microscope  results  agree 
for  agglomerates  in  the  (LW)1/2  between  1  to  2  pm 
range.  The  size  distribution  is  expressed  in  terms  of 
Kg  by  dividing  the  geometric  mean  of  the  (LW)1/2 
for  each  size  class  by  the  factor  2.56.  Measurements 
of  (LW)1/Z  and  R„  for  more  than  285  agglomerates 
produced  a  ratio  between  the  two  parameters  of  2.56 
±  0.33. 

The  number  size  distribution  is  defined  by 


Fn  [log(Kg)] 


dN 

AMlog(Kg) 


AN 

NTAlog(Rg) 


(10) 


where  NT  is  the  total  number  of  analyzed  agglom¬ 
erates,  AN  equal  the  number  of  agglomerates  with 
log(Kg)  between  log(Rgj)  and  log(Kg,:  +  1),  and 
Klog(Kg)  =  log  (Rg  j/Kg  j  +  1).  Fig.  3  provides  a  com¬ 
parison  of  the  number  distribution  for  the  ethene 
and  acetylene  soot  with  the  geometric  mean  Rg  rang¬ 
ing  from  0.02  to  10.9 pm.  It  is  seen  that  the  acetylene 
distribution  is  shifted  to  a  slightly  larger  size  with  a 
geometric  mean  Rg  of  297  nm  versus  217  nm  for  the 
ethene  (see  Table  2).  The  value  of  the  number  dis¬ 
tribution  for  the  largest  size  class  is  about  1000  times 
smaller  than  the  peak  value:  however,  the  large  in¬ 
crease  in  the  total  scattering  cross-section  with  in¬ 
creasing  agglomerate  size  compensates  for  this 
1000-fold  difference  so  that  the  largest  agglomerates 
(even  though  they  are  a  small  fraction  of  the  total 
population)  contribute  significantly  to  the  scattering 
cross-section. 


444 


TURBULENT  COMBUSTION — In-situ  Combustion  Measurements 


TABLE  3 

Uncertainty  analysis  and  correction  factor  for  scattering  measurements  using  TCRN 


Quantity/Effect 

Correction  Factor 

Relative  Standard 

Uncert.  in  pse 

for  543.5  nm  and  632.8  nm  (%) 

Relative  Standard 
Uncert.  in  />se 
for  8.56  nm  (%) 

Linearity  of  PMT 

1.0 

1.5 

Light  extinction  coefficient 

1.7 

2.5 

Attenuation 

0.1 

0.1 

Acceptance  angle 

1.007  (acet.),  1.006  (eth.) 

1.0 

1.0 

Non-ideal  sensor 

1.010 

1.0 

2.0 

Finite  sensor  size 

1.040 

1.0 

1.0 

Total  correction 

1.0.58  (acet.);  1.057  (eth.) 

Total  type  B  uncert.,  UB 

2.6 

3.8 

Total  type  A  uncert..  UA 

1.6 

1.6 

Comb.  rel.  uncert.,  Uc 

3.1 

4.1 

Exp.  rel.  uncert.,  2 Uc 

6.2 

8.2 

Measurement  Uncertainty 

The  measurement  uncertainties  consist  of  type  A 
uncertainties,  which  are  computed  by  statistical 
methods,  and  type  B  uncertainties,  which  involve 
scientific  judgment  rather  than  the  use  of  statistical 
analysis  [12],  Each  uncertainty  component  is  com¬ 
puted  as  relative  standard  uncertainty  defined  as 
the  estimated  standard  deviation  divided  by  the 
mean  value.  For  our  measurements,  the  type  A  un¬ 
certainty  »A,  associated  with  the  mean  value  of  the 
scattering  to  extinction  ratio  pse,  is  computed  based 
on  the  standard  deviation  (SD)  of  the  n  repeat  mea¬ 
surements: 


UA  = 


(11) 


Typically,  five  sets  of  measurements  were  per¬ 
formed  on  each  of  three  days.  The  results  are  in¬ 
cluded  in  Table  3  and  range  between  0.9%  and 
1.6%  of  the  mean  values. 

According  to  equation  8,  the  value  of  pse  is  de¬ 
termined  from  the  photomultiplier  voltage  and  the 
light  extinction  coefficient.  As  discussed  in  Mulhol- 
land  and  Choi  [4],  the  uncertainty  in  the  light  ex¬ 
tinction  coefficient  is  about  1.7%.  For  the  856  nm 
laser,  the  uncertainty  is  greater  and  is  estimated  to 
be  2.5%  due  to  the  relative  lower  response  of  the 
silicon  photodiode  detector  [13],  The  uncertainty  in 
the  linearity  of  the  photomultiplier  output  is  esti¬ 
mated  to  be  1.0%  for  the  visible  wavelengths  and 
1.5%  in  the  near  IR  over  the  typical  range  (factor 
of  three)  in  output  for  the  various  soots  and  dioc- 
tylphthalate  calibration  aerosol. 

Other  factors  affecting  the  type  B  uncertainty  in¬ 
clude  the  limited  acceptance  angle  of  the  detector, 
the  non-ideal  characteristics  of  the  cosine  sensor, 


and  the  finite  size  effect  of  the  detector.  The  influ¬ 
ence  of  limited  acceptance  angle  was  investigated 
by  using  the  size  distribution  (Fig.  3)  and  the 
Fisher-Burford  type  formulation  [3]  for  the  differ¬ 
ential  scattering  cross-section,  Sagg  (Rg,  9): 

Sagg  (Hg>  0)  (12) 

(2nf  d~2Df  Rfl  U  +  cos2  6) 


2 

(4  n 

0\2“ 

1  +  - 

—  fl„ 

sin  - 

3D[ 

U  g 

2/  J 

The  total  scattering,  Sagg(0j,  Of),  integrated  over  the 
angle  range  from  0;  to  Of  and  integrated  over  the 
number  size  distribution  function,  FN[log(Rg)]  from 
log  (Fg.min)  to  log(Rgmax)  is  given  by 

Sjgg  (0i,  Of)  =  f  FN[log(Rg)]d  log(Rg) 

ffff 

^  2 n  sin  0Sagg  (Rg,  0)d0  (13) 

The  ratio  Rscat 


=  sjgg  (ft.  Of) 

Sjgg  (0,  180) 


(14) 


provides  a  measure  of  the  effect  of  the  limited  ac¬ 
ceptance  angle  of  the  detector.  We  find  that  Rscal 
for  the  TCRN  angle  range  of  1.5°  to  178.5°  is  0.993 
for  acetylene  smoke  and  0.994  for  ethene.  The  cor¬ 
responding  correction  factors,  1.007  and  1.006,  are 
used  to  correct  the  measured  scattering  for  acet¬ 
ylene  and  ethene.  The  relative  standard  uncertainty 
in  this  correction  factor  is  estimated  to  be  1.0%. 
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Fig.  4.  Scattering  to  extinction  ratio  for  ethene  and  acet¬ 
ylene  soot  measured  at  543.5  nm,  632.8  nm,  and  856  nm. 
Predictions  were  obtained  using  equation  4.  Previous  mea¬ 
surements  from  Krishnan  et  al.  [7]  in  the  ultraviolet  to  visi¬ 
ble  spectrum  are  also  plotted. 
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Measurements  using  a  goniometer  indicate  that 
the  cosine  sensor/PM  tube  combination  provides 
ideal  response  between  5°  to  175°,  and  that  it  un¬ 
derestimates  the  intensity  at  angles  less  than  5°  or 
greater  than  175°.  The  largest  deviation  of  25%  was 
obtained  at  1°  and  at  179°,  which  are  the  smallest 
and  largest  angles  measured.  These  angular  re¬ 
sponse  characteristics  are  incorporated  in  an  anal¬ 
ysis  similar  to  that  above,  and  the  correction  factor 
and  uncertainty  were  found  to  be  1.010  and  1.0%. 
The  estimates  of  the  effect  of  the  finite  size  of  the 
detector  are  based  on  the  modeling  results  of  Mul- 
holland  and  Bryner  [3],  The  correction  factor  and 
relative  standard  uncertainty  for  finite  detector  size 
are  1.040  and  1%,  respectively. 

The  three  correction  terms  are  independent  of 
agglomerate  sizes  so  that  the  overall  correction 
term  is  obtained  using  the  product  of  the  three 
terms,  with  the  resulting  value  of  1.058  for  acety¬ 
lene  and  1.057  for  ethene.  All  of  the  measured  data 
displayed  in  Table  1  and  Fig.  4  have  been  corrected 
using  these  factors. 

The  results  of  using  the  root  mean  sum  of  squares 
for  combining  the  individual  type  B  uncertainties, 
combining  the  total  type  A  and  type  B  uncertain¬ 
ties,  and  ultimately  computing  the  expanded  rela¬ 
tive  uncertainty,  Uc,  are  summarized  in  Table  3.  Ac¬ 
cording  to  the  definition  of  the  expanded 
uncertainty,  there  is  a  95%  level  of  confidence  that 
pse  is  within  ±  6.2%  of  the  measured  value  for  the 
visible  wavelength  and  ±  8.2%  of  the  measured 
value  for  the  near-IR  wavelengths. 


Discussions 

As  indicated  in  Fig.  4,  the  magnitude  of  the  mea¬ 
sured  pse  (ranging  from  0.196  to  0.316)  is  large 
enough  to  have  a  significant  effect  when  interpret¬ 
ing  light  extinction  data  for  soot  concentration  mea¬ 
surements.  Both  the  measured  and  predicted  values 
of  pse  (predictions  using  equation  4)  at  543.5  nm, 
632.8  nm,  and  856  nm  are  plotted  in  Fig.  4.  The 
input  values  for  the  predictions  (using  equation  4) 
are  taken  from  measurements  of  soot  structure 
from  Table  2  and  a  refractive  index  value  of  1.55 
+  i0.8  [2],  We  focus  on  the  trends  regarding  the 
effects  of  primary  sphere  size  and  wavelength 
rather  than  the  quantitative  agreement  because  of 
the  uncertainty  in  the  refractive  index  of  the  soot. 
For  both  measurement  and  prediction,  the  acety¬ 
lene  soot  (which  has  a  larger  average  primary  par¬ 
ticle  size)  has  the  larger  value  of  pse  compared  to 
ethene.  This  same  effect  is  apparent  in  Fig.  4  for 
the  results  of  Krishnan  et  al.  [7],  The  predicted  ra¬ 
tio  between  acetylene  and  ethene  pse,  approxi¬ 
mately  1.6,  at  632.8  nm  is  larger  than  the  measured 
ratio  of  approximately  1.15.  This  discrepancy  may 
be  due  to  beam  shielding  effects  in  which  attenu¬ 
ation  of  light  by  the  primary  particles  on  the  front 
side  of  the  agglomerate  reduces  the  intensity  reach¬ 
ing  the  rest  of  agglomerate  [13], 

As  shown  in  Fig.  4,  the  measured  change  in  pse 
with  wavelength  is  in  qualitative  agreement  with  the 
predictions  for  543.5  and  632.8  nm.  However,  there 
is  a  large  difference  at  856  nm,  where  the  data  in¬ 
dicate  almost  no  change  in  pse  at  856  nm  compared 
to  633  nm,  while  the  predicted  value  decreased  by 
about  35%  for  both  ethene  and  acetylene  soot.  This 
invariance  in  the  scattering  ratio  is  an  important 
finding,  because  it  is  often  thought  that  the  scatter¬ 
ing  in  the  near-IR  becomes  less  important  than  in 
the  visible  spectrum  (see  equation  4).  These  are  the 
first  quantitative  measurements  of  scattering  to  ex¬ 
tinction  ratio  in  the  IR.  One  potential  explanation 
is  a  wavelength  dependence  of  the  soot  refractive 
index.  Although  additional  experiments  at  longer 
wavelengths  will  be  required  to  clarify  this  matter, 
the  present  measurements  have  provided  important 
information  regarding  the  scattering  behavior  in  the 
visible  and  near-IR  spectrum. 

Fig.  4  also  displays  pse  measurements  at  351  nm, 
400  nm,  488  nm,  514  nm,  and  632.8  nm  reported 
by  Krishnan  et  al.  [7]  for  overfire  soot  from  a  5  cm 
turbulent  acetylene  and  ethene  flames.  The  mea¬ 
sured  values  of  pse  in  the  visible  spectrum  by  Krish¬ 
nan  et  al.  [7]  increases  with  wavelength,  while  our 
data  indicate  a  decrease  over  the  visible  wavelength 
range.  The  pse  measurements  of  Krishnan  et  al.  [7] 
at  514.5  nm  are  within  10%  of  the  present  mea¬ 
surements  at  543.5  nm.  However,  measurements  at 
632.8  nm  for  both  acetylene  and  ethene  for  Krish¬ 
nan  et  al.  [7]  are  more  than  50%  larger  than  the 
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present  measurements.  The  estimated  expanded 
uncertainty  (95%  confidence  level)  of  ±6.2%  for 
the  present  data  is  much  smaller  than  the  differ¬ 
ence  in  the  results  for  the  two  studies.  Wu  et  al. 

[6] ,  who  used  the  same  method  as  Krishnan  et  al. 

[7] ,  report  expanded  relative  uncertainties  of  10% 
for  total  scattering  and  5%  for  total  extinction,  lead¬ 
ing  to  an  expanded  uncertainty  of  the  ratio  pse  equal 
to  about  11%.  Also,  the  differences  between  the 
two  sets  of  experiments  are  not  entirely  due  to  the 
different  burner  configurations  (laminar  burner  in 
the  present  study  versus  turbulent  burner  used  in 
Krishnan  et  al.  [7])  since  comparisons  of  the  present 
laminar  flame  measurements  with  measurements  by 
Mulholland  and  Choi  [4]  obtained  using  the  TCRN 
and  the  same  turbulent  flame  of  Krishnan  et  al.  [7] 
were  within  15%.  The  acceptance  angle  range  for 
the  scattering  measurements  by  Krishnan  et  al.  [7] 
is  5°  to  160°.  The  effect  of  this  limited  angle  range 
based  on  our  measured  size  distribution  is  a  reduc¬ 
tion  of  the  measured  scattering  by  about  8%.  This 
has  the  effect  of  increasing  the  difference  between 
the  present  measurements  and  those  of  Krishnan  et 
al.  [7],  While  more  study  is  required  to  resolve  the 
differences  in  the  measured  pse,  both  studies  show 
that  scattering  component  is  a  significant  portion  of 
the  extinction  cross-section. 
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